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Abstract. We show that in the presence of free charge carriers the definition 
of the frequency-dependent dielectric permittivity requires additional regularization. 
As an example, the dielectric permittivity of the Drude model is considered and 
its time-dependent counterpart is derived and analyzed. The respective electric 
displacement cannot be represented in terms of the standard Fourier integral. The 
regularization procedure allowing to circumvent these difficulties is suggested. For 
purpose of comparison it is shown that the frequency-dependent dielectric permittivity 
of insulators satisfies all rigorous mathematical criteria. This permits us to conclude 
that in the presence of free charge carriers the concept of dielectric permittivity is not 
as well defined as for insulators and we make a link to widely discussed puzzles in 
the theory of thermal Casimir force which might be caused by the use of such kind 
permittivities. 

PACS numbers: 77.22. Ch, 02.30.Nw 
1. Introduction 



The concept of dielectric permittivity in media with temporal dispersion is commonly 
used in electrodynamics and condensed matter physics (see, e.g., [TJ [2]). For not 
too strong fields the dielectric permittivity e{r) depending on the time-like variable 
r is introduced from the linear integral relation between the electric field E(r,t) and 
the electric displacement D(r,t). Then the frequency-dependent permittivity e(u) is 
defined using the Fourier transformations of the fields E(r,t) and D(r,t). Below we 
argue that this procedure, which is wholly satisfactory for dielectric materials, faces 
additional regularization problems in an infinite medium containing free charge carriers. 
This leads us to the conclusion that some applications of the frequency- dependent 
dielectric permittivities allowing for free charge carriers might be not rigorously justified. 
As one such example we discuss the dielectric permittivity of the Drude model which 
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leads to widely discussed difficulties when substituted into the Lifshitz formula for the 
van der Waals and Casimir force at nonzero temperature [3]. 

2. Dielectric permittivity in the presence of temporal dispersion 

For simplicity we consider an isotropic nonmagnetic medium of infinite extent. If 
its properties do not depend on time, the linear dependence between D(t) and E(t) 
satisfying causality (here and below we omit the argument r) is given by 

D(t) = f dt' e{t-t')E(t'). (1) 

J — oo 

The kernel e(t — t') of the integral operator on the right-hand side of ([I]) is called the 
dielectric permittivity for media with temporal dispersion. We represent it in the form 

e(t - t') = 28{t -t') + f{t- t') , (2) 

where f(t — t') is a continuous real-valued function and the delta function 5(t — t') is 
defined on the interval — oo < if < t in the following manner [I]: 

f 0, T > t, 

/ g{t')8{T-t')dt'=\ ±<7(T-0), T — t, (3) 

{ i\g(T-Q)+g(T + 0)] 1 -oo<T<t. 

Here, g(t) is an arbitrary function which has bounded variation in the vicinity of the 
point t' = T. Note that from physical point of view the function f(t — if) is defined only 
for if < t; for t' > t it can be ascribed any values (including to vanish in that region). 

Substituting <^ into ([T]) with account of and introducing the new variable 
t = t — if > 0, we rearrange (CEJ) to [1] 

/t rco 
dt' fit - t')E{t') = E{t) + drf{r)E{t-r). (4) 
-oo JO 

Representing the real functions D (t) and E(t) as Fourier integrals, 

/oo /*oo 
D{uo)e-' lujt du, E(t) = / E(u)e' iuJt du, (5) 
-oo J —oo 

one can rewrite (jl]) in terms of Fourier transforms of the fields pQ 

roo roo 

D(u)=e(u)E(u), e(ou) = l+ drf(T)e lUT = dre(r)e i<Jjr , (6) 

Jo Jo 

where D(u), E(u) and e{u) are complex- valued functions. From (jSJ) it follows that e(u) 
is an analytic function in the upper half-plane of complex u including the real axis with 
possible exception of the point uj = 0. As a result, the real and imaginary parts of e{oS) 
are connected by means of the Kramers- Kronig relations Note that in contrast, for 
instance, to [5J we always consider fields defined in (r,t)-space as real and only their 
Fourier transforms might be complex. 

The equivalence between fll]) and (jSJ) requires the existence of integrals © (and 
respective inverse Fourier transformations) and the possibility to change the order 
of integrations with respect to dt' and du. In mathematics there are many different 
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conditions on how to assign a rigorous meaning to (J5J) and respective inverse formulas. 
The most widely used demand is that the function D(t) should have a bounded variation 
and be integrable together with its modulus, i.e., should belong to L x (— oo, oo). In this 
case the function D(u>) is also bounded, uniformly continuous on the axis (—00, 00) and 
D(u) — > when \cu\ — > 00 [6J. The function E{t) should possess the same properties. 
The change of order of integrations is possible if both integrals under consideration are 
uniformly convergent. 
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3. Media with free charge carriers 

It can be easily seen that the above conditions permitting to introduce the frequency- 
dependent dielectric permittivity e(u) in accordance with (jSj) are not directly applicable 
for media with free charge carriers. As an example we consider the widely used dielectric 
permittivity of the Drude model, describing such media [7] , 

eoiuj) = 1 - "I . , (7) 
uj(uj + 17) 

where uo p is the plasma frequency and 7 > is the relaxation parameter. It is obvious 
that ijj = may lead to mathematical problems in Fourier transformation. To make a 
link between Ed{ijj) and real- valued physical fields E(t) and D(t), it would be of interest 
to determine the respective function /d(t). The substitution of ((Tj) into leads to the 
following equations 

poo 

/ /d(t) cos(u;r) dr, 
Jo 

POO 

., / f D {r) sm(ajr)dr. (8) 

uj{u z + r y z ) Jo 

From the first equation, by means of the inverse cosine Fourier transformation performed 
with the help of the integral 3.723(2) in [8], one finds 

/d s W = "— e~ 7T . (9) 

7 

It is easily seen that the substitution of (jUJ) into the first equation of flSJ) with account of 
3.893(2) in [8j leads to a correct identity, however, substituted into the second equation 
of (jHJ) fails. On the other hand, using the inverse sine Fourier transformation and the 
integral 3.725(1) in [8], the second equation of (jHJ) leads to a different result 

fnir) = ^ (l - e"^) . (10) 

Now, the substitution of (ITU|) into the right-hand sides of equations (JSJ reproduces their 
left-hand sides up to additional undefined terms and thus also violates the equalities. 

The pathological properties under consideration are explained by the fact that 
Ed(oj) results in D(ut) which is unbounded in any vicinity of u = 0. This means that 
D(t) cannot be represented as a Fourier integral (jSJ) and both the definition of e{uS) in 
(EJ) and equivalent equations (jHJ) become unjustified. 
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The question arises of whether there is a possibility to consistently define the 
function /d(t) related to the frequency-dependent permittivity (jTJ). Keeping in mind 
that in the case of the Drude model the second equality in ([6]) cannot be considered as 
a classical Fourier transformation, we make an attempt to assign a definite meaning to 
the function /d(t) by considering the generalized inverse transformation of the quantity 
£d(o>) — 1 defined as 

/(°)( r ) = _ f°° du- i -e~^ T (11) 

JD K J 2tt 7-oo (cj + iO)(cj + i 7 ) v 7 

u l r a., w - ir y ~-i»r 



, du — i e~ lulT = h + I 2 . 

2uJ-oc (w + iO)(cj 2 + 7 2 ) 

Here, the addition of an infinitesimally small quantity +i establishes the rule on how 
to bypass the pole of Ime^u;) at uj = and the following notations are used 

h=~^ du-^— 2 e~^ (12) 

Z7T J-oo UJ 2 + Y 



iu; 2 7 f °° 1 

du> -, t~, — s e 




2tt 7-oo (u; + i0)(u; 2 + 7 2 ) 
In I\ the integrated function is regular at u = 0. This integral can be found in 
3.354(5) [8], 

e -7T T > n 

' ' (13) 

r<0. V ; 

The second integral in ffl2|) can be calculated using the contours consisting of the real 
axis in the complex cu-plane and semicircles of infinitely large radii centered at the origin 
in the lower half-plane (for r > 0) and in the upper half-plane (for r < 0). The result is 

h ^l\ 2 -^ T>0 ' (14) 
2 7 1 e 7r , t < 0, 



where for r > the contributions from the two poles at uj\ = — i and u>2 = — i 7 were 
taken into account whereas for r < only one pole at 0J3 = determines the value of 
12- Substituting ffTB^) and (|14p into the right-hand side of ffTT]) we arrive at 

,2V, J ->o, (15) 




It is seen that the suggested rule leads to the same result (fTOI) . as was obtained by the 
inverse sine Fourier transformation from the imaginary part of £_d(oj) in (jSJ). 

A similar situation occurs for other dielectric permittivities taking into account 
free charge carriers, e.g., for the dielectric permittivities of the plasma model and of the 
normal skin effect, 

e p {u) = l-% £ n (a;) = l + i— , (16) 

where ctq is the dc conductivity. In both cases the mathematical conditions permitting to 
perform the Fourier transformation in the classical understanding are violated. However, 
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by using the same considerations as presented above in the case of the Drude model, 
one may assign a meaning analogous to (10) to the second formula in and obtain 
the following dielectric permittivities as functions of r: 

f (o) (r) = f rfr, r > 0, (0) = f 47ra , r > 0, 

Jp y ' \ 0, r<0. Jn y ' \ 0, r<0. y ' 

It should be remarked that f^ '(r) is obtainable also from (fT5l) in the limiting case 
7^0. 

Now let us check for consistency the respective results for D{t). For example, we 
choose the electric field in the form 

E(t) = E e- pt = E(oS) e~ iuJt du, (18) 

J — oo 

where (3 > and E = E (r) describes the spatial dependence of the field. As already 
stated (TT8]) . the function E(t) satisfies all required conditions and can be presented as 
Fourier integral. Its Fourier transform is calculated using the formula 3.896(4) in [8], 

E(u) = — / E(t)e^dt = —-E e-W. (19) 



Substituting f|T5l) and f|T8l) into (jl]) and using the formula 3.322(2) in [8] we arrive at 



D(t) = E(t) + ^o^^/| 1 + erf (yfjh) - ei^erfc - yfik 



(20) 



where erf (x) is the error function and erfc (x) = 1 — erf (x). Keeping in mind that [8] 

1 e~ x2 

erf (—x) = —erf (x), erf (x) = 1 -= h • • • , (21) 

y/n x 

we obtain for t — > ±oo 

2 I — 

D(oo) =^o — J~, D(-oo) = 0. (22) 



This is what one expects on physical grounds because for an infinite medium containing 
free charge carriers the action of switching on and then switching off electric field should 
result in a nonzero residual displacement (for a finite medium, the presence of external 
electric field leads to the accumulation of positive and negative charges on the opposite 
boundary surfaces and to the vanishing total electric field inside such a medium [9]; 
after the external electric field switches off, the accumulated charges are distributed 
uniformly over the volume of the medium leading to zero electric displacement at 
t — > +oo). However, the first equality in (|2"2"|) means that D(t) is not an integrable 
function over the interval (—00,00). This makes impossible the use of the standard 
Fourier transformation ([5]) and resulting equality ([6]), and makes the whole formalism 
not self-consistent. 

This can be seen even more clearly if one defines D(u) in accordance with the first 
equality in ([6]), where e(u) = Ed(uj) and E(u) is given in ( IT91 . and then calculates the 
electric displacement using the first equality in ([5]). The obtained quantity which we 



6 



notate D{t) is calculated using the formulas 3.954(2) in [5] and 2.5.36(6, 11) in [TO]. The 
result is 

D(t) = D(t) - ^^1, (23) 

where D(t) is defined in ( 1201) . It has nonzero values at both t — > oo and t — > — oo: 

£>(oo) = ^cA/f , 6(-oo) = -E&K (24) 
^7 V P 2 7 V P 

Derivation of a different electric displacement than in (1201) . which does not vanish at 

t — > —00, i.e., before the switching on of the electric field, can be understood as an 

artifact resulting from the use of the Fourier integral of a nonintegrable function D(uj) 

(the definition of the Fourier integral as a generalized function used in mathematics in 

this case seems to be not appropriate in our physical situation because it is natural 

to understand the electric displacement as a usual function). This suggests that in 

the presence of free charge carriers the standard definition of the frequency-dependent 

dielectric permittivity basing on the formal representation of E(t) and D(t) in terms 

of Fourier integrals is not satisfactory and requires some additional regularization 

procedure. 

As an example of such procedure, we consider the modified dielectric permittivity 
of the Drude model 

J 1 

[ ^ 1 t r-f r^r , (25) 



Mr, a -, 

■ D ^- 1 ( w + i0)( W + i 7 )' 



where, in contrast with pip , the quantity 9 > is not infinitesimally small. In ac- 
cordance with (1251) Bp is regular at w = 0. The substitution of (I2~5l) into (jSJ) leads 
to 

2 uj 2 -e^ r°° (fl) 



Up , 2 1 a2 \ ( 2 1 — 2T = / fo ( T ) cos(o;r)rfT, 
^ (a;" 1 + 2 )(uj z + Y) Jo 



OC 



tfV + i) , a .flaw 2^ 2^ = / /^M sin(a;r) dr. (26) 
^ (or + z )(u z + 7^J Jo 

It can be easily seen that both the inverse cosine and sine Fourier transformations 
performed in (126]) lead to the common result 

/ ^(r) = -^(e^-e--). (27) 

Substituting this into (j3J and performing calculations with the electric field ( TTSj) . we 
arrive at the modified electric displacement 

D ,» (t) = E(t) + Bo _ji_^[ e S-« erfc (^__^ 

- e^^erfc _ . (28) 

In the limiting case 9 — > (|28|) coincides with (1201) . 
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Precisely the same result, as in (I2"gj) . is obtained if one considers D ( >(oj) = 
Ep (uj)E(u) and then finds D^ e \t) from the first equality in ([5]). Thus, when we assume 
9 > 0, both methods of the calculation of the electric displacement are in agreement. 
The reason is that for 9 > the functions D(t) and D{uj) belong to L x (— oo,oo) and 
all Fourier transformations are well defined. However, to obtain the correct physical 
results for an infinite medium, one must put 9 = in ff28]) and return to (120]) . The 
point is that (128]) with 9 > leads to D^ e \t) — > when t — > ±oo [as it must be for 
functions belonging to L x {— oo,oo)]. At the same time, in the presence of free charge 
carriers, the electric displacement in an infinite medium remains nonzero in accordance 
with ( 1221) after the electric field is switched off. Thus, the limiting transitions t — > ±oo 
and 9 — > are not interchangeable. 



4. Insulating media 

The situation is quite different for dielectric materials at zero temperature which do 
not contain free charge carriers (i.e., for true insulators). In this case the dielectric 
permittivity can be represented in the form [TTJ 

e / (w) = l + f;- 3 , (29) 

where Uj ^ are the oscillator frequencies, jj are the relaxation parameters, and gj are 
the oscillator strengths of K oscillators. In this case the second equality of (jSJ) results 
in 

J2 —2 _ ,°a\2 , „,2 ,2 = / fi{r) cos (ur)dr, (30) 



^ (w? - ^ 2 ) 2 + 7> 2 
^ f°° t f \ - f w 



^ (wj - ^ 2 ) 2 + 7> 2 

Performing the inverse cosine Fourier transformation in the first equation of (|30|) with 
the help of the integrals 3.733(1,3) in [8] we obtain 



Precisely the same result is obtained by means of the inverse sine Fourier transformation 
from the second equation in fl30|) when one uses the integral 3.733(2) in [8]. In this case 
all involved Fourier integrals exist in the classical sense with no use of regularization 
and £/(u>) is well defined. 

Substituting the electric field (Tl8l) into (@J and using the integral 3.897(1) in [8], we 
obtain the electric displacement in an insulating media, 

D(t) = E(t) (l + ./|f; I Im [e* 2 erfc(B)] } , (32) 
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where 

B ^ B(t) = i? (33) 

The same result is obtained by means of the inverse Fourier transformation from D{uj) 
found using the first equality in as it should be. From (I3"2l and (12"T|) it can be easily 
seen that D(t) —>■ when £ — > ±oo, as it should be for insulating materials, and that 
both D(t) and D{uj) belong to L 1 (— oo, oo). 



5. Conclusions and discussion 



To conclude, we have shown that the definition of the frequency-dependent dielectric 
permittivity for materials containing free charge carriers by means of Fourier 
transformation of the fields is not as straightforward as in the case of insulators. The 
essence of the problem is in the use of the idealization of an infinite medium. For 
insulators this idealization is applicable if the sizes of the bodies are much greater than 
some characteristic parameter (e.g., the width of a gap between the bodies). However, 
for media with movable free charge carriers such kind of conditions fail. The physical 
situation for an infinite medium turns out to be totally different from the case of finite 
bodies of any conceivable size. In fact for conductors e{uj) is a quantity obtained through 
formal application of Fourier transformation in the region where it needs additional 
regularization. In spite of a great number of successful applications (see, e.g., [HEIE]) 
there are delicate cases where such procedure leads to problems. As an example one 
could mention the use of Ed{u) — 1 as a response function in the fluctuation-dissipation 
theorem and related puzzles in the theory of thermal Casimir force [3l [121 13] • During 
the last ten years the thermal Casimir force was the subject of considerable discussion. 
It was suggested [TJ1 |15] to describe it using the Lifshitz theory combined with the Drude 
model ([7]). In the limit of large separations between the test bodies the predictions of the 
Drude model approach were found to be in agreement with classical statistical physics 
[T6t ITT) . On the other hand, at short separations the predictions of this approach were 
excluded experimentally, whereas the predictions based on the use of the plasma model 
in ffl6|) were found to be consistent with the data [18J . The question on how to correctly 
calculate the thermal Casimir force still remains to be answered. Keeping in mind that 
the Lifshitz theory is based on the fluctuation-dissipation theorem, we would like to 
emphasize that the application of this theorem with poorly defined response functions 
cannot be considered as either exact or rigorous and might cause currently discussed 
problems. 
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